IMAGE CODING METHOD AND APPARATUS 
AND IMAGE DECODING METHOD AND APPARATUS 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to an image processing 
technology, and more particularly relates to an image coding 
and decoding method and apparatus therefor which takes error 
processing into consideration. 

2. Description of the Related Art 

Recently, image processing and compression methods such 
as those proposed by MPEG (Motion Picture Expert Group) have 
expanded to be used with transmission media such as network 
and broadcast rather than just storage media such as CDs. 
Generally speaking, the success of the digitization of 
broadcast materials has been caused at least in part by the 
availability of MPEG compression coding technology. In this 
way, a barrier that previously existed between broadcast and 
other types of communication has begun to disappear, leading 
to a diversification of service-providing businesses. Thus, 
we are facing a situation where it is hard to predict how the 
digital culture would evolve in this age of broadband. 
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Even in such a chaotic situation, it is clear that the 
direction of the compression technology of motion pictures 
will be to move to both higher compression rates and better 
image quality. It is a well-known fact that block distortion 
5 in MPEG compression is sometimes responsible for causing 
degraded image quality and preventing the compression rate 
from being improved. Another problem that occurs in 
compression algorithms is poor tolerance of errors. Unlike 
reading data from storage media, this is particularly an issue 
10 for transmission of data through networks such as the Internet 
in which data delivery can be subject to a greater number of 
errors . 

15 SUMMARY OF THE INVENTION 

The present invention has been made in view of the 
foregoing circumstances and an object thereof is to provide a 
new compression technology for motion pictures which also 
20 provides error tolerance. 

In the following, an image decoding technology and an 
image coding technology according to the present invention are 
described primarily with regard to motion pictures. However, 
the use of the technology is not limited to motion pictures 



and includes other image processing, for example, image 
effects such as morphing and walkthrough. 

A preferred embodiment according to the present 
invention relates to an image decoding apparatus. First, for 
5 nonnegative integers i and j , when data of an i-th key frame 
and a j-th key frame are defined as KFi and KFj, respectively, 
a correspondence data file between the i-th key frame and the 
j-th key frame is defined as Ci,j, and data of an intermediate 
frame between the i-th key frame and the j-th key frame is 

10 defined as IFi,j. Now, the image decoding apparatus includes: 
an error detector which receives a data stream that includes 
KFi, KFi+i and Ci,i+i and which detects whether or nor there is 
an error in the data stream; an intermediate image generator 
which generates IFi,i+i from the data stream; and an error 

15 controller which, when an error occurs in the data stream, 
controls said intermediate image generator in a manner such 
that an error avoidance processing is performed in the 
intermediate image generator. Hereinafter, the term "frame", 
which is a unit of image, will not be distinguished from the 

20 term "image" unless otherwise necessary. Further, the term 

"error" indicates a data error and other arbitrary faults and 
obstruction. 

The intermediate image generator may generate iFi^j by 
interpolating using Ci,j, KFi and KFj. IFi,j may correspond to 



an arbitrary frame between KFi and KFj, and, in that sense, 
IFi,j representatively symbolizes all intermediate frames that 
may be generated between KFi and KFj . 

Another preferred embodiment according to the present 
invention relates to an image decoding method. This method 
includes: acquiring a data stream that includes a plurality of 
key frames and a correspondence data file therebetween; 
generating an intermediate frame between the key frames from 
the data stream; and monitoring for an error in the data 
stream, wherein, when an error is detected, an error avoidance 
processing is performed at the time of generating the 
intermediate frame. 

In the above apparatus and method, there may be several 
types of error avoidance processing, depending on where the 
error occurs, for example, as follows: 

1. The intermediate frame is generated from data of another 
key frame which is substituted for data of the error- 
containing key frame. 

2. The intermediate frame is generated based on data of 
another key frame and a correspondence data file relating 
thereto. 

3. The intermediate frame is generated from another 
correspondence data file which is substituted for the error- 
containing correspondence data file. 



4. The intermediate frame is generated based on a new 
correspondence data file generated based on another 
correspondence data file. 

Moreover, when it is judged that a seriousness of the 
error is below a predetermined level, the error avoidance 
processing may be abandoned. Situations in which the serious 
of the error is low may include the following examples: 

1. In all pixels of a key frame, the total number of error- 
containing pixels is low. 

2. In the image of a key frame, the error is not present or 
substantially less in the low frequency component thereof. 

3. Even though there are a plurality of error-containing 

pixels or areas, they are relatively dispersed. 

4. The error-containing area is close to the edge of an 
image . 

5. The error-containing area is in the background of an 
object and so forth, thus being inconspicuous. 

6. The intensity of an error-containing pixel or area is low. 

Still another preferred embodiment according to the 
present invention relates to an image coding apparatus. The 
image coding apparatus includes: an image input unit which 
receives data of key frames; a correspondence data generator 
which generates Ci,i+i by utilizing KFi and KFi+i, and generates 
Ci,j by utilizing KFi and KFj, where j > i+1, among the input 



key frame data ; and a stream generator which generates a data 
stream including data generated by the correspondence data 
generator . 

Still another preferred embodiment according to the 
present invention relates to an image coding method. The 
image coding method includes: first generating, based on data 
of two adjacent key frames, correspondence data therebetween; 
second generating, based on data of key frames disposed at 
some interval, correspondence data therebetween; and 
generating a data stream in a manner such that the 
correspondence data generated in said first generating serves 
as main data and the correspondence data generated in said 
second generating serves as spare data. 

Here, the term "main data" indicates data to be 
reproduced in a usual manner while the term "spare data" 
indicates data to be used for repairing an error which has 
occurred at the time of decoding processing. Thus, the second 
generating may be executed with lower frequency than execution 
of the first generating. Moreover, in the present method, the 
correspondence data may be generated in a manner such that 
data for checking for errors is embedded therein. An example 
of a type of error checking data is a parity bit which is 
inserted into the correspondence data. 

The image coding method according to the present 



invention may further include, when generating the 
correspondence data, an image matching which includes: 
multiresolutionalizing two key frames by respectively 
extracting critical points; performing, for example, a pixel- 
by-pixel matching computation on those key frames between 
similar multiresolution levels; and acquiring a pixel-by-pixel 
correspondence relation at a finest level of resolution while 
inheriting a result of the pixel-by-pixel matching computation 
at a different multiresolution level. 

In the above-described embodiments, the matching method 
that utilizes critical points may be an application of the 
technology (referred to as the "base technology" hereinafter) 
proposed by Japanese Patent No. 2927350 which is owned by the 
same assignee as the present patent specification. 

It is to be noted that it is also possible to have 
replacement or substitution of the above-described structural 
components and elements of methods in part or whole as between 
method and apparatus or to add elements to either method or 
apparatus and also, the apparatuses and methods may be 
implemented by a computer program and saved on a recording 
medium or the like and are all effective as and encompassed by 
the present invention. 

Moreover, this summary of the invention includes 
features that may not be necessary features such that an 



embodiment of the present invention may also be a sub- 
combination of these described features. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1(a) is an image obtained as a result of the 
application of an averaging filter to a human facial image. 

Fig. 1(b) is an image obtained as a result of the 
application of an averaging filter to another human facial 
image. 

Fig. 1(c) is an image of a human face at p*^'°' obtained 
in a preferred embodiment in the base technology. 

Fig. 1(d) is another image of a human face at p'^'°' 
obtained in a preferred embodiment in the base technology. 

Fig. 1(e) is an image of a human face at p*^'"'-' obtained 
in a preferred embodiment in the base technology. 

Fig. 1(f) is another image of a human face at p'^'-"-' 
obtained in a preferred embodiment in the base technology. 

Fig. 1(g) is an image of a human face at p*^'^' obtained 
in a preferred embodiment in the base technology. 

Fig. 1(h) is another image of a human face at p'^'^' 
obtained in a preferred embodiment in the base technology. 

Fig. l(i) is an image of a human face at p'^'"^^ obtained 



in a preferred embodiment in the base technology. 

Fig. l(j) is another image of a human face at p^^'^' 
obtained in a preferred embodiment in the base technology. 
Fig. 2 (R) shows an original quadrilateral. 
5 Fig. 2(A) shows an inherited quadrilateral. 

Fig. 2(B) shows an inherited quadrilateral. 
Fig. 2(C) shows an inherited quadrilateral. 
!^ Fig. 2(D) shows an inherited quadrilateral. 

^ Fig. 2(E) shows an inherited quadrilateral. 

10 Fig. 3 is a diagram showing the relationship between a 

l^Z source image and a destination image and that between the m-th 

level and the (m-l)th level, using a quadrilateral. 

U 

Fig. 4 shows the relationship between a parameter?? 
(represented by x-axis) and energy Cf (represented by y-axis) . 
15 Fig. 5(a) is a diagram illustrating determination of 

whether or not the mapping for a certain point satisfies the 
bijectivity condition through the outer product computation. 

Fig. 5(b) is a diagram illustrating determination of 
whether or not the mapping for a certain point satisfies the 
20 bijectivity condition through the outer product computation. 

Fig. 6 is a flowchart of the entire procedure of a 
preferred embodiment in the base technology. 

Fig. 7 is a flowchart showing the details of the process 
at SI in Fig. 6. 
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Fig. 8 is a flowchart showing the details of the process 
at SIO in Fig. 7. 

Fig. 9 is a diagram showing correspondence between 
partial images of the m-th and (m-l)th levels of resolution. 
5 Fig. 10 is a diagram showing source hierarchical images 

generated in the embodiment in the base technology. 

Fig. 11 is a flowchart of a preparation procedure for S2 
in Fig. 6. 

Fig. 12 is a flowchart showing the details of the 
10 process at S2 in Fig. 6. 

Fig. 13 is a diagram showing the way a submapping is 
determined at the 0-th level. 

Fig. 14 is a diagram showing the way a submapping is 
determined at the first level. 
15 Fig. 15 is a flowchart showing the details of the 

process at S21 in Fig. 12. 

Fig. 16 is a graph showing the behavior of energy 
corresponding to f^""'^^ (A. = iAA) which has been obtained for a 
certain f t""'^' while varying A . 
20 Fig. 17 is a diagram showing the behavior of energy Cy"^ 

corresponding to f'"' (7 = /A//) (r = 0,1,...) which has been obtained 
while varying rj . 

Fig. 18 shows a structure of an image coding apparatus 
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according to an embodiment of the invention. 

Fig, 19 shows a relationship between a main file and a 
spare file generated by a matching processor. 

Fig. 20 shows a structure of a data stream generated by 
a stream generator. 

Fig. 21 shows another relationship between a main file 
and a spare file. 

Fig. 22 shows another structure of a data stream 
generated by the stream generator. 

Fig. 23 shows still another relationship between a main 
file and a spare file. 

Fig. 24 shows still another structure of a data stream 
generated by the stream generator. 

Fig. 25 shows an area of a key frame which may be used 
for checking for any errors during decoding. 

Fig. 2 6 shows an area of a correspondence data file 
which may be used for checking for any errors during decoding. 

Fig. 27 shows a structure of an image decoding apparatus 
according to an embodiment of the invention. 

Fig. 28 shows a principle of generation of intermediate 
frames by an intermediate image generator. 

Fig. 2 9 shows an error avoidance processing performed 
when an error is detected in a key frame. 

Fig. 30 shows another error avoidance processing 
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performed when an error is detected in a key frame. 

Fig, 31 shows an error avoidance processing performed 
when errors are detected in key frames. 

Fig. 32 shows still another error avoidance processing 
5 performed when an error is detected in a key frame. 

Fig. 33 shows an error avoidance processing performed 
when an error is detected in a correspondence data file. 

Fig. 34 shows another error avoidance processing 
performed when an error is detected in a correspondence data 



DETAILED DESCRIPTION OF THE INVENTION 

15 The invention will now be described based on the 

preferred embodiments, which do not intend to limit the scope 
of the present invention, but exemplify the invention. All of 
the features and the combinations thereof described in the 
embodiment are not necessarily essential to the invention. 

20 First, the multiresolutional critical point filter 

technology and the image matching processing using the 
technology, both of which will be utilized in the preferred 
embodiments, will be described in detail as "Base Technology". 
Namely, the following sections [1] and [2] (below) belong to 



the base technology, where section [1] describes elemental 
techniques and section [2] describes a processing procedure. 
These techniques are patented under Japanese Patent No. 
2927350 and owned by the same assignee of the present 
invention. However, it is to be noted that the image matching 
techniques provided in the present embodiments are not limited 
to the same levels. In particular. In Figs. 18 to 34, image 
data coding and decoding techniques utilizing, in part, the 
base technology will be described in more detail. 

Base Technology 

[1] Detailed description of elemental techniques 

[1.1] Introduction 

Using a set of new multiresolutional filters called 
critical point filters, image matching is accurately computed. 
There is no need for any prior knowledge concerning the 
content of the images or objects in question. The matching of 
the images is computed at each resolution while proceeding 
through the resolution hierarchy. The resolution hierarchy 
proceeds from a coarse level to a fine level. Parameters 
necessary for the computation are set completely automatically 
by dynamical computation analogous to human visual systems. 
Thus, There is no need to manually specify the correspondence 
of points between the images . 



The base technology can be applied to, for instance, 
completely automated morphing, object recognition, stereo 
photograrametry, volume rendering, and smooth generation of 
motion images from a small number of frames. When applied to 
morphing, given images can be automatically transformed. When 
applied to volume rendering, intermediate images between cross 
sections can be accurately reconstructed, even when a distance 
between cross sections is rather large and the cross sections 
vary widely in shape. 

[1.2] The hierarchy of the critical point filters 

The multiresolutional filters according to the base 
technology preserve the intensity and location of each 
critical point included in the images while reducing the 
resolution. Initially, let the width of an image to be 
examined be N and the height of the image be M. For 
simplicity, assume that N=M=2n where n is a positive integer. 
An interval [0, N] c R is denoted by I . A pixel of the image 
at position (i, j) is denoted by p*""'^* where i,j E I. 

Here, a multiresolutional hierarchy is introduced. 
Hierarchized image groups are produced by a multiresolutional 
filter. The multiresolutional filter carries out a two 
dimensional search on an original image and detects critical 
points therefrom. The multiresolutinal filter then extracts 
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the critical points from the original image to construct 
another image having a lower resolution. Here, the size of 
each of the respective images of the m-th level is denoted as 
2"'X2"' (0 ^ m ^ n) . A critical point filter constructs the 
following four new hierarchical images recursively, in the 
direction descending from n. 

/,;«f =max(min(pS5^ 
;7(\«f =min(max(pS5^^ 

(1) 

where we let 

— (2) 

The above four images are referred to as subimages 
hereinafter. When itiinxst^x+i and maxx^t^x+i are abbreviated to ot 
and /S , respectively, the subimages can be expressed as 
15 follows: 

P^"--'^ =a{x)P{y)p^"'^'''^ 

Namely, they can be considered analogous to the tensor 
products of a and )S . The subimages correspond to the 
respective critical points. As is apparent from the above 
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equations, the critical point filter detects a critical point 
of the original image for every block consisting of 2 X 2 
pixels. In this detection, a point having a maximum pixel 
value and a point having a minimum pixel value are searched 
5 with respect to two directions, namely, vertical and 
horizontal directions, in each block. Although pixel 
intensity is used as a pixel value in this base technology, 
various other values relating to the image may be used. A 
pixel having the maximum pixel values for the two directions, 

10 one having minimum pixel values for the two directions, and 
one having a minimum pixel value for one direction and a 
maximum pixel value for the other direction are detected as a 
local maximum point, a local minimum point, and a saddle 
point, respectively. 

15 By using the critical point filter, an image (1 pixel 

here) of a critical point detected inside each of the 
respective blocks serves to represent its block image (4 
pixels here) in the next lower resolution level. Thus, the 
resolution of the image is reduced. From a singularity 

20 theoretical point of view, a (x) a (y) preserves the local 

minimum point {minima point) , (x) fi (y) preserves the local 
maximum point (maxima point) , a {k) 0 (y) and/S (x) a (y) preserve 
the saddle points. 

At the beginning, a critical point filtering process is 



applied separately to a source image and a destination image 
which are to be matching-computed. Thus, a series of image 
groups, namely, source hierarchical images and destination 
hierarchical images are generated. Four source hierarchical 

5 images and four destination hierarchical images are generated 
corresponding to the types of the critical points. 

Thereafter, the source hierarchical images and the 
destination hierarchical images are matched in a series of 
resolution levels. First, the minima points are matched using 

10 p'"''"^ Next, the first saddle points are matched using p^"''^^ 
based on the previous matching result for the minima points. 
The second saddle points are matched using p'"''^^. Finally, 
the maxima points are matched using p'""'^'. 

Figs. Ic and Id show the subimages p*^'°^ of the images in 

15 Figs, la and lb, respectively. Similarly, Figs, le and If 
show the subimages p^^'^', Figs. Ig and Ih show the subimages 
p'^'^>, and Figs, li and Ij show the subimages p<^'^' . 
Characteristic parts in the images can be easily matched using 
subimages. The eyes can be matched by p'^'°' since the eyes 

20 are the minima points of pixel intensity in a face. The 
mouths can be matched by p'^'^' since the mouths have low 
intensity in the horizontal direction. Vertical lines on both 
sides of the necks become clear by p'^'^' . The ears and bright 
parts of the cheeks become clear by p^^'^^ since these are the 
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maxima points of pixel intensity. 

As described above, the characteristics of an image can 
be extracted by the critical point filter. Thus, by 
comparing, for example, the characteristics of an image shot 
by a camera with the characteristics of several objects 
recorded in advance, an object shot by the camera can be 
identified. 

[1.3] Computation of mapping between images 

Now, for matching images, a pixel of the source image at 
the location (i,j) is denoted by pf")) and that of the 
destination image at (k,l) is denoted by where i, j, k, 1 

e I. The energy of the mapping between the images (described 
later in more detail) is then defined. This energy is 
determined by the difference in the intensity of the pixel of 
the source image and its corresponding pixel of the 
destination image and the smoothness of the mapping. First, 
the mapping f :p«-'°> ^ q(-'°' between p'^'°> and q^"''^^ with the 
minimum energy is computed. Based on f^"''^', the mapping f^^^'^^ 
between p*""'^' and q'""'^* with the minimum energy is computed. 
This process continues until f'"''^' between p'"*'^' and q^"''^^ is 
computed. Each f'™'^' (i = 0,1,2,...) is referred to as a 
submapping. The order of i will be rearranged as shown in the 
following equation (3) in computing f'"*'^* for reasons to be 
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. described later. 

j-irrtJ) . pimM')) _^ qimM')) ( 3 ) 

where cr (i) e {0,1,2,3}. 

5 [1. 3. 1] Bijectivity 

When the matching between a source image and a 
destination image is expressed by means of a mapping, that 
f% mapping shall satisfy the Bijectivity Conditions (BC) between 

ifl the two images (note that a one-to-one surjective mapping is 

~'~ 10 called a bijection) . This is because the respective images 

should be connected satisfying both surjection and injection, 
'»3 and there is no conceptual supremacy existing between these 

H images. It is to be noted that the mappings to be constructed 

-° here are the digital version of the bijection. In the base 

15 technology, a pixel is specified by a co-ordinate point. 

The mapping of the source sub image (a sub image of a 
source image) to the destination subimage (a subimage of a 
destination image) is represented by f*'"'^': 1/2''""' X 1/2"""' 
I/a^-^'X 1/2"-"' (s = 0,1,...), where /(^f =(^,0 means that ^i^f of 
20 the source image is mapped to q'^^^^ of the destination image. 
For simplicity, when f (i,j) = (k,l) holds, a pixel q(k,i) is 
denoted by qf (i,j) - 

When the data sets are discrete as image pixels (grid 



points) treated in the base technology, the definition of 
bij activity is important. Here, the bijection will be defined 
in the following manner, where i, j, k and 1 are all integers. 
First, a square region R defined on the source image plane is 

considered 

(•n,s) (m,s) (m,s) (m,s) 

P(J,j) P{i+l,j)P(,+\j+ViP{i,m ^ ^ ' 

where i = 0, 2"'-l, and j = 0, 2"-l. The edges of R are 
directed as follows: 



{in,s)„(m,s) „(m,i) „(/n^) „(.m,s) (m,s) j (/»,*) {m,s) ,c\ 



This square region R will be mapped by f to a 

quadrilateral on the destination image plane: 

This mapping f <'*'2' (R) , that is, 

should satisfy the following bijectivity conditions (referred 
to as BC hereinafter) : 

1. The edges of the quadrilateral f (R) should not 

intersect one another. 

2. The orientation of the edges of f (R) should be the same 
as that of R (clockwise in the case shown in Fig. 2, described 
below) . 

3. As a relaxed condition, a retraction mapping is allowed. 

Without a certain type of a relaxed condition as in, for 
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example, condition 3 above, there would be no mappings which 
completely satisfy the BC other than a trivial identity 
mapping. Here, the length of a single edge of f (R) may be 
zero. Namely, f '""'^^ (R) may be a triangle. However,^ f (R) 

5 is not allowed to be a point or a line segment having area 
zero. Specifically speaking, if Fig. 2R is the original 
quadrilateral. Figs. 2A and 2D satisfy the BC while Figs 2B, 
2C and 2E do not satisfy the BC. 

In actual implementation, the following condition may be 

10 further imposed to easily guarantee that the mapping is 

surjective- Namely, each pixel on the boundary of the source 
image is mapped to the pixel that occupies the same location 
at the destination image. In other words, f (i,j)=(i,j) (on 
the four lines of i=0, i=2"'-l, j=0, j=2"'-l) . This condition 

15 will be hereinafter referred to as an additional condition. 

[1. 3. 2] Energy of mapping 

[1. 3. 2. 1] Cost related to the pixel intensity 

The energy of the mapping f is defined. An objective 
20 here is to search a mapping whose energy becomes minimum. The 
energy is determined mainly by the difference in the intensity 
between the pixel of the source image and its corresponding 
pixel of the destination image. Namely, the energy C^^fy^ of 
the mapping f at (i, j ) is determined by the following 
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equation (7) . 

where ) and V(q%fj)) are the intensity values of the 

pixels /jjjf and , respectively. The total energy C*"''^' of 

f is a matching evaluation equation, and can be defined as the 
sum of Cj";f as shown in the following equation (8) . 

/=0 j=0 

[1. 3. 2. 2] Cost related to the locations of the pixel for 
smooth mapping 

In order to obtain smooth mappings, another energy Df 
for the mapping is introduced. The energy Df is determined by 
the locations of and (i=0, 1, 2'^-l, j=0, 1, 2"^-!) , 

regardless of the intensity of the pixels. The energy £){";)^ of 
the mapping f*"''^' at a point (i,j) is determined by the 
following equation (9). 

where the coefficient parameter 77 which is equal to or 
greater than 0 is a real number. And we have 

^oS"?=iio''i)-/'"'^^o-,7)ir — (10) 



<'=/-u'=y-i 

where 

\\(x,y)\\ = ^x'+y' (12), 

i' and j ' are integers and f(i',j') is defined to be zero for 
i'<0 and j''<0. Eo is determined by the distance between {i,j) 
and f (i, j) . Eo prevents a pixel from being mapped to a pixel 
too far away from it. However, as explained below, Eo can be 
replaced by another energy function. Ei ensures the 
smoothness of the mapping. Ei represents a distance between 
the displacement of p(i,j) and the displacement of its 
neighboring points. Based on the above consideration, another 
evaluation equation for evaluating the matching, or the energy 
Df is determined by the following equation: 

[1. 3. 2. 3] Total energy of the mapping 

The total energy of the mapping, that is, a combined 
evaluation equation which relates to the combination of a 
plurality of evaluations, is defined as MJ^f-'^ + Df'''\ where A 
> 0 is a real number. The goal is to detect a state in which 
the combined evaluation equation has an extreme value, namely, 
to find a mapping which gives the minimum energy expressed by 
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the following: 

rmaiXCf-'^ + Df'"^} — (14) 

Care must be exercised in that the mapping becomes an 
identity mapping if A=0 and ?? =0 (i.e., f '"''^Mi, j ) = (i/ j ) for 
5 all i=0, l^-.^a""-! and j=0, 1, 2"'-l) . As will be described 

later, the mapping can be gradually modified or transformed 
from an identity mapping since the case of A =0 and r? =0 is 
evaluated at the outset in the base technology. If the 
combined evaluation equation is defined as Cj"'*^ + IDj"'*^ where 

10 the original position of A is changed as such, the equation 
with A=0 and 77=0 will be Cj"'*^ only. As a result thereof, 
J pixels would randomly matched to each other only because their 

pixel intensities are close, thus making the mapping totally 
meaningless. Transforming the mapping based on such a 

15 meaningless mapping makes no sense. Thus, the coefficient 
parameter is so determined that the identity mapping is 
initially selected for the evaluation as the best mapping. 

Similar to this base technology, differences in the 
pixel intensity and smoothness are considered in a technique 

20 called "optical flow" that is known in the art. However, the 
optical flow technique cannot be used for image transformation 
since the optical flow technique takes into account only the 
local movement of an object. However, global correspondence 
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can also be detected by utilizing the critical point filter 
according to the base technology. 



[1. 3. 3] Determining the mapping with multiresolution 



5 



A mapping fnd.n which gives the minimum energy and 



satisfies the BC is searched by using the multiresolution 
hierarchy. The mapping between the source sub image and the 
destination subimage at each level of the resolution is 
computed. Starting from the top of the resolution hierarchy 

10 (i.e., the coarsest level), the mapping is determined at each 
resolution level, and where possible, mappings at other levels 
are considered. The number of candidate mappings at each 
level is restricted by using the mappings at an upper (i.e., 
coarser) level of the hierarchy. More specifically speaking, 

15 in the course of determining a mapping at a certain level, the 
mapping obtained at the coarser level by one is imposed as a 
sort of constraint condition. 

We thus define a parent and child relationship between 
resolution levels. When the following equation (15) holds. 
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where |_xj denotes the largest integer not exceeding x. 



are respectively called the parents of _P(^"y 



and ^{"jf,. Conversely, p^f and are the child of p'^^^;}^ 



and the child of 4"!/)''^' respectively. A function parent (i,j) 
is defined by the following equation (16) : 

pare«^(/,j) = ([^J|^J) (16) 

Now, a mapping between ^j? determined by 

computing the energy and finding the minimum thereof. The 
value of f ^""'^Mi, j) = (k/l) is determined as follows using f (m- 
l^s) (m=l, 2,...,n) . First of all, a condition is imposed that 
^Jj)^ should lie inside a quadrilateral defined by the 
following definitions (17) and (18) . Then, the applicable 
mappings are narrowed down by selecting ones that are thought 
to be reasonable or natural among them satisfying the BC. 

(m,s) (m,s) {m,s) (y) ^7) 

where 

gC".-) y) = iparentii, j)) + (jjarent(i, j) + (1,1)) — (18) 

The quadrilateral defined above is hereinafter referred 
to as the inherited quadrilateral of p^^jf . The pixel 
minimizing the energy is sought and obtained inside the 
inherited quadrilateral. 

Fig. 3 illustrates the above-described procedures. The 
pixels A, B, C and D of the source image are mapped to A' , B' , 
C and D' of the destination image, respectively, at the (m- 
l)th level in the hierarchy. The pixel p^^J^^ should be mapped 
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to the pixel which exists inside the inherited 

quadrilateral A'B'C'D'. Thereby, bridging from the mapping at 
the (m-l)th level to the mapping at the ra-th level is 
achieved. 

5 The energy Eq defined above may now be replaced by the 

following equations (19) and (20) : 

^oov)=||/^'"'°H/,7)-g^'"HU-)|r — (19) 

^oa,) ^l/'^'CUV/'^'C^^f'Ci^O — (20) 

for computing the submapping f*"*'"' and the submapping f*""'^' at 

10 the m-th level, respectively. 

In this manner, a mapping which maintains a low energy 
of all the submappings is obtained. Using the equation (20) 
makes the submappings corresponding to the different critical 
points associated to each other within the same level in order 

15 that the subimages can have high similarity. The equation 
(19) represents the distance between f^"''^'(i,j) and the 
location where (i,j) should be mapped when regarded as a part 
of a pixel at the (m-l)the level. 

When there is no pixel satisfying the BC inside the 

20 inherited quadrilateral A'B'C'D', the following steps are 
taken. First, pixels whose distance from the boundary of 
A'B'C'D' is L (at first, L=l) are examined. If a pixel whose 
energy is the minimum among them satisfies the BC, then this 



pixel will be selected as a value of f (i, j ) . L is 
increased until such a pixel is found or L reaches its upper 
bound L^^. is fixed for each level m. If no pixel is 

found at all, the third condition of the BC is ignored 
temporarily and such mappings that caused the area of the 
transformed quadrilateral to become zero (a point or a line) 
will be permitted so as to determine f '"''^^ (i, j ) • If such a 
pixel is still not found, then the first and the second 
conditions of the BC will be removed. 

Multiresolution approximation is essential to 
determining the global correspondence of the images while 
preventing the mapping from being affected by small details of 
the images. Without the multiresolution approximation, it is 
impossible to detect a correspondence between pixels whose 
distances are large. In the case where the multiresolution 
approximation is not available, the size of an image will 
generally be limited to a very small size, and only tiny 
changes in the images can be handled. Moreover, imposing 
smoothness on the mapping usually makes it difficult to find 
the correspondence of such pixels. That is because the energy 
of the mapping from one pixel to another pixel which is far 
therefrom is high. On the other hand, the multiresolution 
approximation enables finding the approximate correspondence 
of such pixels. This is because the distance between the 



pixels is small at the upper (coarser) level of the hierarchy 
of the resolution. 

[1. 4] Automatic determination of the optimal parameter values 

One of the main deficiencies of the existing image 
matching techniques lies in the difficulty of parameter 
adjustment. In most cases, the parameter adjustment is 
performed manually and it is extremely difficult to select the 
optimal value. However, according to the base technology, the 
optimal parameter values can be obtained completely 
automatically. 

The systems according to this base technology include 
two parameters, namely, A and ?? , where A and ?? represent 
the weight of the difference of the pixel intensity and the 
stiffness of the mapping, respectively. In order to 
automatically determine these parameters, the are initially 
set to 0. First, A is gradually increased from A =0 while rj 
is fixed at 0 . As A becomes larger and the value of the 
combined evaluation equation (equation (14) ) is minimized, the 
value of Cj^^ for each submapping generally becomes smaller. 
This basically means that the two images are matched better. 
However, if A exceeds the optimal value, the following 
phenomena occur: 

1. Pixels which should not be corresponded are erroneously 
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corresponded only because their intensities are close. 

2. As a result, correspondence between images becomes 
inaccurate, and the mapping becomes invalid. 

3. As a result, Df''^ in equation (14) tends to increase 
5 abruptly. 

4. As a result, since the value of equation (14) tends to 
increase abruptly, f^'^'^^ changes in order to suppress the 

i abrupt increase of Df-'K As a result, Cf'"^ increases. 

Therefore, a threshold value at which Cf-"^ turns to an 

I y 

!%■> 

■Q 10 increase from a decrease is detected while a state in which 

111 

- equation (14) takes the minimum value with A being increased 

o 

is kept. Such A is determined as the optimal value at 77=0. 

Next, the behavior of Cf-"^ is examined while r? is increased 

gradually, and r? will be automatically determined by a method 
15 described later. A will then again be determined 

corresponding to such an automatically determined 77 . 

The above-described method resembles the focusing 

mechanism of human visual systems. In the human visual 

systems, the images of the respective right eye and left eye 
20 are matched while moving one eye. When the objects are 

clearly recognized, the moving eye is fixed. 



[1. 4. 1] Dynamic determination of A 



Initially, A is increased from 0 at a certain interval, 
and a subimage is evaluated each time the value of A changes. 
As shown in equation (14), the total energy is defined by 
^^(m,*) _^ ^ equation (9) represents the smoothness 

5 and theoretically becomes minimum when it is the identity 
. mapping. Eq and Ei increase as the mapping is further 
distorted. Since Ei is an integer, 1 is the smallest step of 
Df-"^. Thus, it is impossible to change the mapping to reduce 
^ the total energy unless a changed amount (reduction amount) of 

10 the current ACj"f is equal to or greater than 1. Since Df''^ 
'I ■ increases by more than 1 accompanied by the change of the 

mapping, the total energy is not reduced unless ^C^Tj) is 
J; reduced by more than 1 . 

Under this condition, it is shown that C["f^ decreases in 
15 normal cases as A increases. The histogram of C^"}^ is 

denoted as h(l), where h(l) is the number of pixels whose 
energy C^"f is 1^. In order that A 1^ ^ 1 for example, the 
case of 1^=1/ A is considered. When A varies from Ai to X2r 
a number of pixels (denoted A) expressed by the following 
20 equation (21) : 



— (21) 
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changes to a more stable state having the energy shown in 
equation (22) : 

Cf''^ -f= Cf''^ -J. — ( 22 ) 

Here, it is assumed that the energy of these pixels is 
5 approximated to be zero. This means that the value of C^^j^^ 
changes by: 

_ QCO»^)=-A (23) 

''E As a result, equation (24) holds, 

i 3Cr_m — (24) 

10 Since h(l)>0 , Cj"'*^ decreases in the normal case. However, 
C3 when A exceeds the optimal value, the above phenomenon, that 

.St is, an increase in C^f'"^ occurs. The optimal value of A is 

determined by detecting this phenomenon. 
When 

15 h{l) = Hl'^-^ — (25) 

is assumed, where both H{H>0) and k are constants, the 
equation (26) holds: 

dX ~ X"'^"^ ^^^^ 
Then, if k5^-3, the following equation (27) holds: 



The equation (27) is a general equation of Cj™'*^ (where C is a 
constant) . 

When detecting the optimal value of A , the number of 
pixels violating the BC may be examined for safety. In the 
course of determining a mapping for each pixel, the 
probability of violating the BC is assxomed as a value po here. 
In this case, since 

^^'^ 

holds, the nxamber of pixels violating the BC increases at a 

rate of: 

Thus, 

^^ = 1 (30) 

PoKD 

is a constant. If it is assumed that h(l)=Hl'^, the following 
equation (31) , for example, 

B,A"'''"=p,H — (31) 
becomes a constant. However, when A exceeds the optimal 
value, the above value of equation (31) increases abruptly. 
By detecting this phenomenon, i.e. whether or not the value of 
Bf)^'^^'''^ /2"' exceeds an abnormal value B^^^^^, the optimal value 



of A can be determined. Similarly, whether or not the value 
of B^?^'^*^'"^ 12"" exceeds an abnormal value can be used to 

check for an increasing rate Bi of pixels violating the third 
condition of the BC. The reason why the factor 2"* is 
5 introduced here will be described at a later stage. This 

system is not sensitive to the two threshold values -S^^^^^^ and 
^xthres ' threshold values ^o/Are^ ^^'^ ^Uhres used to 

detect excessive distortion of the mapping which may not be 
v^i detected through observation of the energy Cy™'*^ . 

)3 10 In the experimentation, when A exceeded 0.1 the 

^ computation of f*"*'^' was stopped and the computation of f'^'^+i' 

was started. That is because the computation of submappings 
is affected by a difference of only 3 out of 255 levels in 
pixel intensity when A > 0.1 and it is then difficult to 
15 obtain a correct result. 



[1. 4. 2] Histogram h(l) 

The examination of Cj™'*^ does not depend on the histogram 
h(l), however, the examination of the BC and its third 
20 condition may be affected by h(l). When (A , Cj"'*M is 

actually plotted, k is usually close to 1. In the experiment, 
k=l is used, that is, B^?^ and are examined. If the true 

value of k is less than 1, B^^ and B^A^ are not constants and 



increase gradually by a factor of Z^^'""^'^ . If h(l) is a 
constant, the factor is, for example, J^''^ . However, such a 
difference can be absorbed by setting the threshold Bq^,^^^ 
appropriately . 

Let us model the source image by a circular object, with 
its center at(xo,yo) and its radius r, given by: 



with its center at(xi,yi) and radius r. In the above, let 
c(x) have the form of c(x)=x''. When the centers (xo,yo) and 
(xi,yi) are sufficiently far from each other, the histogram 
h(l) is then in the form: 

h(l)ozrl''(k^O) — (34) 

When k=l, the images represent objects with clear 
boundaries embedded in the background. These objects become 
darker toward their centers and brighter toward their 
boundaries. When k=-l, the images represent objects with 
vague boundaries. These objects are brightest at their 




— (32) 



and the destination image given by: 




— (33) 
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centers, and become darker toward their boundaries. Without 
much loss of generality, it suffices to state that objects in 
images are generally between these two types of objects. 
Thus, choosing k such that -l^k^l can cover most cases and the 
5 equation (27) is generally a decreasing function for this 
range . 

As can be observed from the above equation (34), 
attention must be directed to the fact that r is influenced by 
'r the resolution of the image, that is, r is proportional to 2". 

10 This is the reason for the factor 2" being introduced in the 
'I^ above section [1.4.1]. 

[1. 4. 3] Dynamic determination of 77 

The parameter tj can also be automatically determined in 

15 a similar manner. Initially, rj is set to zero, and the final 
mapping f'"^ and the energy C|"^ at the finest resolution are 
computed. Then, after r) is increased by a certain value Atj , 
the final mapping f''^' and the energy Cj"^ at the finest 
resolution are again computed. This process is repeated until 

20 the optimal value of 77 is obtained. 77 represents the 
stiffness of the mapping because it is a weight of the 
following equation (35) : 

K£=\\f''"''\U)-f""''-'\i,j)f — (35) 
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If 7? is zero, Dj"-* is determined irrespective of the 
previous subinapping, and the present submapping may be 
elastically deformed and become too distorted. On the other 
hand, if 77 is a very large value, D^f is almost completely 
5 determined by the immediately previous submapping. The 

submappings are then very stiff, and the pixels are mapped to 
almost the same locations. The resulting mapping is therefore 
the identity mapping. When the value of 77 increases from 0, 
Cj"^ gradually decreases as will be described later. However, 

10 when the value of 7? exceeds the optimal value, the energy 

starts increasing as shown in Fig. 4. In Fig. 4, the x-axis 
represents 77, and y-axis represents Cf. 

The optimum value of 7? which minimizes Cf can be 
obtained in this manner. However, since various elements 

15 affect this computation as compared to the case of A , 

changes while slightly fluctuating. This difference is caused 
because a submapping is re-computed once in the case of A 
whenever an input changes slightly, whereas all the 
submappings must be re-computed in the case of r? . Thus, 

20 whether the obtained value of Cf is the minimum or not cannot 
be determined as easily. When candidates for the minimum 
value are found, the true minimum needs to be searched by 
setting up further finer intervals. 
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[1. 5] Supersampling 

When deciding the correspondence between the pixels, the 
range of f'"^'^' can be expanded to R X R (R being the set of 
5. real numbers) in order to increase the degree of freedom. In 
this case, the intensity of the pixels of the destination 
image is interpolated, to provide f*™'^' having an intensity at 
non-integer points: 

= 10 That is, supersampling is performed. In an example 

implementation, f'"''^' may take integer and half integer 
Q values, and 

'J ^(?£?(o.s.o.5)) — <37) 

^ is given by 

15 iViqlrjf) + V(qi-fi,,,,))/2 (38) 



[ 1. 6] Normalization of the pixel intensity of each image 

When the source and destination images contain quite 
different objects, the raw pixel intensity may not be used to 
20 compute the mapping because a large difference in the pixel 

intensity causes excessively large energy Cj™'*'-' and thus making 
it difficult to obtain an accurate evaluation. 

For example, a matching between a human face and a cat's 



face is computed as shown in Figs. 20(a) and 20(b). The cat's 
face is covered with hair and is a mixture of very bright 
pixels and very dark pixels. In this case, in order to 
compute the submappings of the two faces, subimages are 
5 normalized. That is, the darkest pixel intensity is set to 0 
while the brightest pixel intensity is set to 255, and other 
pixel intensity values are obtained using linear 
interpolation. 

10 [1. 7] Implementation 

In an example implementation, a heuristic method is 
utilized wherein the computation proceeds linearly as the 
source image is scanned. First, the value of f*""'^^ is 
determined at the top leftmost pixel (i,j)=(0,0). The value 

15 of each f'"''®'(i,j) is then determined while i is increased by 
one at each step. When i reaches the width of the image, j is 
increased by one and i is reset to zero. Thereafter, 
f"^'^'(i,j) is determined while scanning the source image. 
Once pixel correspondence is determined for all the points, it 

20 means that a single mapping f*™-^' is determined. 

When a corresponding point qf(i,j) is determined for P(i,j), 
a corresponding point qf(i,j+i) of P(i,j+i) is determined next. 
The position of qf(i,j+i) is constrained by the position of qf(i,j) 
since the position of qf{i,j+i) satisfies the BC. Thus, in this 
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system, a point whose corresponding point is determined 
earlier is given higher priority. If the situation continues 
in which (0,0) is always given the highest priority, the final 

mapping might be unnecessarily biased. In order to avoid this 
5 bias, f'™'=' is determined in the following manner in the base 
technology. 

First, when (s mod 4) is 0, f'"*'^' is determined starting 
from (0,0) while gradually increasing both i and j. When (s 
mod 4) is 1, f*"^'^' is determined starting from the top 

10 rightmost location while decreasing i and increasing j . When 
(s mod 4) is 2, f'""'^' is determined starting from the bottom 
rightmost location while decreasing both i and j . When (s mod 
4) is 3, f''"'^' is determined starting from the bottom leftmost 
location while increasing i and decreasing j . Since a concept 

15 such as the submapping, that is, a parameter s, does not exist 
in the finest n-th level, f*™'^* is computed continuously in two 
directions on the assumption that s=0 and s=2 . 

In this implementation, the values of f'"^'^*(i,j) 
(m=0,...,n) that satisfy the BC are chosen as much as possible 

20 from the candidates (k,l) by imposing a penalty on the 

candidates violating the BC. The energy D(k,i) of a candidate 
that violates the third condition of the BC is multiplied by 
0 and that of a candidate that violates the first or second 
condition of the BC is multiplied by ^ . In this 



implementation, <t)=2 and ^=100000 are used. 

In order to check the above-mentioned BC, the following 
test may be performed as the procedure when determining 
(k, 1) =f (i, j ) . Namely, for each grid point (k,l) in the 
inherited quadrilateral of f*™'^*(i,j), whether or not the z- 
component of the outer product of 

W = AxB (39) 

is equal to or greater than 0 is examined, where 

^ = —(40) 

B = '^^A~^ —(41) 
Here, the vectors are regarded as 3D vectors and the z-axis is 
defined in the orthogonal right-hand coordinate system. When 
W is negative, the candidate is imposed with a penalty by 
multiplying D^'^jf by t/j so that it is not as likely to be 
selected. 

Figs. 5(a) and 5(b) illustrate the reason why this 
condition is inspected. Fig. 5(a) shows a candidate without a 
penalty and Fig. 5(b) shows one with a penalty. When 
determining the mapping f ^"^'^^ (i, j+1) for the adjacent pixel at 
(i,j+l), there is no pixel on the source image plane that 
satisfies the BC if the z-component of W is negative because 
then ^J;,^ passes the boundary of the adjacent quadrilateral. 
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[1. 7. 1] The order of submappings 

In this implementation, o-(0)=0, o-(l)=l, cr {2)=2, 
cr(3)=3, o-(4)=0 are used when the resolution level is even, 
while o-(0)=3, cr(l)=2, ct(2)=1, o-{3)=b, a-(4)=3 are used when 
5 the resolution level is odd. Thus, the submappings are 
shuffled to some extent. It is to be noted that the 
submappings are primarily of four types, and s may be any of 0 
to 3. However, a processing with s=4 is used in this 

■'Z implementation for a reason to be described later. 

'& 10 

'.^ [1. 8] Interpolations 

After the mapping between the source and destination 
images is determined, the intensity values of the 
corresponding pixels are interpolated. In the implementation, 

X 15 trilinear interpolation is used. Suppose that a square 

P(i,j)P(i+i,j>P(i+i, j+i)P(i, j+i) on the source image plane is mapped to 
a quadrilateral qf (i, j)qf (i+i,j)qf (i+i, j+i)qf {i,j+i) on the destination 
image plane. For simplicity, the distance between the image 
planes is assumed to be 1. The intermediate image pixels 
20 r(x,y,t) (0:^x<N-l, 0<y^M-l) whose distance from the source 

image plane is t (0:^t^l) are obtained as follows. First, the 
location of the pixel r(x,y,t), where x,y,teR, is determined 
by equation (42): 



(X, y) = (\- dx)(l - dy){\ - t){i, j) + (1 - dx)i\ - dy)tf{i, j) 

+ dx{\-dy\\-t)(i + \,j) + dx{\-dy)tf{i + \J) ___ 
+ (1 - dx)dy{\ - t)(i, / + 1) + (1 - dx)dytf(i, 7+1) 
+ dxdy(\ - tXi + 1, 7 + 1) + dxdytfii + + 1) 

The value of the pixel intensity at r(x,y,t) is then 

determined by equation (43) : 

V(r{x,yj)) = (1 - dx)i\ - dy){\ - t)V{p,,^j,) + (1 - dx){\ - dyWiq^u)) 
+ dx(l - dy)(l - t)V(p,,,,j, ) + dx(l- dy)tV{q,,,,,^^, ) 
+ (1 - dx)dy(\ - 0F(/7(,,,i)) + (1 - dx)dytViqf^,^^,,^) 
+ dxdyil - Of^(/'(,.i.,.i) ) + dxdytV{qf^,,,j,,) ) 

(43) 

where dx and dy are parameters varying from 0 to 1 . 

[1. 9] Mapping to which constraints are imposed 

So far, the determination of a mapping in which no 
constraints are imposed has been described. However, if a 
correspondence between particular pixels of the source and 
destination images is provided in a predetermined manner, the 
mapping can be determined using such correspondence as a 
constraint . 

The basic idea is that the source image is roughly 
deformed by an approximate mapping which maps the specified 
pixels of the source image to the specified pixels of the 
destination image and thereafter a mapping f is accurately 
computed . 

First, the specified pixels of the source image are 
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mapped to the specified pixels of the destination image, then 
the approximate mapping that maps other pixels of the source 
image to appropriate locations are determined. In other words, 
the mapping is such that pixels in the vicinity of a specified 
5 pixel are mapped to locations near the position to which the 
specified one is mapped. Here, the approximate mapping at the 
m-th level in the resolution hierarchy is denoted by f'™' . 

The approximate mapping F is determined in the following 
manner. First, the mappings for several pixels are specified. 
10 When Us pixels 

PihJ^\p(}xJi\"->Pi\-iJn..-i) — (44) 

of the source image are specified, the following values in the 
equation (45) are determined. 

F^''\i„j\) = (kM,.-, (45) 

15 For the remaining pixels of the source image, the amount 

of displacement is the weighted average of the displacement of 
p(ih^jh) (h=0,..., ng-l) . Namely, a pixel P(i,j) is mapped to the 
following pixel (expressed by the equation (46) ) of the 
destination image. 

J^(k„-i^,l^-j^)weight^ii,J) 
20 F^'"\i,j) = ^ — — (4 6) 

where 



weight^(i,j) = -f- . , ".. (47) 

total _weight(i,j) 

where 

total _weight(hj) = ''Y!m(h-hj\-j)f (48) 

h=0 

Second, the energy D^"^^ of the candidate mapping f is 
changed so that a mapping f similar to F^"*' has a lower energy. 
Precisely speaking, ^(,"')^ is expressed by the equation (4 9) : 



\If^"'^ (i, j) - Z^™'*^ (i, j)f , otherwise 



— (50) 



where K,p>0. Finally, the resulting mapping f is determined 
by the above-described automatic computing process. 

Note that E^"^^^ becomes 0 if f (i, j ) is sufficiently 
close to F^""(ifj) i.e., the distance therebetween is equal to 
or less than 



,2(«-.)J 



(51) 



This has been defined in this way because it is desirable to 
determine each value f'"''^*(i,j) automatically to fit in an 

appropriate place in the destination image as long as each 
value f'"''^'(i,j) is close to F*"''(i,j). For this reason, there 
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is no need to specify the precise correspondence in detail to 
have the source image automatically mapped so that the source 
image matches the destination image. 

5 [2] Concrete Processing Procedure 

The flow of a process utilizing the respective elemental 
techniques described in [1] will now be described. 

Fig. 6 is a flowchart of the overall procedure of the 
base technology. Referring to Fig. 6, a source image and 

10 destination image are first processed using a 

multiresolutional critical point filter (SI) . The source image 
and the destination image are then matched (S2) . As will be 
understood, the matching (S2) is not required in every case, 
and other processing such as image recognition may be 

15 performed instead, based on the characteristics of the source 
image obtained at SI. 

Fig. 7 is a flowchart showing details of the process Si 
shown in Fig. 6. This process is performed on the assumption 
that a source image and a destination image are matched at S2. 

20 Thus, a source image is first hierarchized using a critical 
point filter (SIO) so as to obtain a series of source 
hierarchical images. Then, a destination image is hierarchized 
in the similar manner (Sll) so as to obtain a series of 
destination hierarchical images. The order of SIO and Sll in 
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the flow is arbitrary, and the source image and the 
destination image can be generated in parallel. It may also be 
possible to process a number of source and destination images 
as required by subsequent processes. 
5 Fig. 8 is a flowchart showing details of the process at 

SIO shown in Fig. 7. Suppose that the size of the original 
source image is 2^X2^. Since source hierarchical images are 
sequentially generated from an image with a finer resolution 
to one with a coarser resolution, the parameter m which 

10 indicates the level of resolution to be processed is set to n 
(SlOO) . Then, critical points are detected from the images 
p(m,o)^ p<™'i\ p(m,2) pti^'^^of the m-th level of resolution, 

using a critical point filter (SlOl) , so that the images p^"'" 
1,0)^ p(m-i,i)^ p(m-i,2) p(m-i,3) (m-l)th level are 

15 generated (8102). Since m=n here, p*'"'^' ^p*"''^' =pC»'2) =p(m.3) 

=p*"' holds and four types of subimages are thus generated from 
a single source image. 

Fig. 9 shows correspondence between partial images of 
the m-th and those of (m-l)th levels of resolution. Referring 

20 to Fig. 9, respective number ic values shown in the figure 
represent the intensity of respective pixels. p'^irs) 
symbolizes any one of four images p*"''^' through p'™'"^', and when 
generating pC^-^'O'^ p(m,o) used from p*""'^' . For example, as 
for the block shown in Fig. 9, comprising four pixels with 
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their pixel intensity values indicated inside, images p'-^-^'°^ ^ 
p(m-i,i)^ p(m-i,2) p(m-i,3) g^quire -3", -8", - 6" and "10", 

respectively, according to the rules described in [1.2]. This 
block at the m-th level is replaced at the (m-l)th level by 
respective single pixels thus acquired. Therefore, the size 
of the subimages at the (m-l)th level is 2"'"^X2"'"^. 

After m is decremented {S103 in Fig. 8), it is ensured 
that m is not negative (S104) . Thereafter, the process 
returns to SlOl, so that subimages of the next level of 
resolution, i.e., a next coarser level, are generated. The 
above process is repeated until subimages at m=0 (0-th level) 
are generated to complete the process at SIO. The size of the 
subimages at the 0-th level is 1 X 1. 

Fig. 10 shows source hierarchical images generated at 
SIO in the case of n=3. The initial source image is the only 
image common to the four series followed. The four types of 
subimages are generated independently, depending on the type 
of critical point. Note that the process in Fig. 8 is common 
to Sll shown in Fig. 7, and that destination hierarchical 
images are generated through a similar procedure. Then, the 
process at SI in Fig. 6 is completed. 

In this base technology, in order to proceed to S2 shown 
in Fig. 6 a matching evaluation is prepared. Fig. 11 shows 
the preparation procedure. Referring to Fig. 11, a plurality 



of evaluation equations are set (S30) . The evaluation 
equations may include the energy Cf''^ concerning a pixel 
value, introduced in [1.3.2.1], and the energy D^"'*^ concerning 
the smoothness of the mapping introduced in [1.3.2.2]. Next, 
by combining these evaluation equations, a combined evaluation 
equation is set (S31) . Such a combined evaluation equation 
may be XC^J^jf + Df'"'' . Using r? introduced in [1.3.2.2], we have 

YL(^\;:f^^vE:^::i^E^£) — (52) 

In the equation (52) the sum is taken for each i and j where i 
and j run through 0, 1,... , 2"* Now, the preparation for 
matching evaluation is completed. 

Fig. 12 is a flowchart showing the details of the 
process of S2 shown in Fig. 6. As described in [1], the 
source hierarchical images and destination hierarchical images 
are matched between images having the same level of 
resolution. In order to detect global correspondence 
correctly, a matching is calculated in sequence from a coarse 
level to a fine level of resolution. Since the source and 
destination hierarchical images are generated using the 
critical point filter, the location and intensity of critical 
points are stored clearly even at a coarse level. Thus, the 
result of the global matching is superior to conventional 
methods . 
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Referring to Fig. 12, a coefficient parameter rj and a 
level parameter m are set to 0 (S20) . Then, a matching is 
computed between the four subimages at the m-th level of the 
source hierarchical images and those of the destination 
5 hierarchical images at the m-th level, so that four types of 
submappings f^""'^^ (s=0, 1, 2, 3) which satisfy the BC and 
minimize the energy are obtained (S21) . The BC is checked by 
■using the inherited quadrilateral described in [1.3.3]. In 
that case, the submappings at the m-th level are constrained 

10 by those at the (m-l)th level, as indicated by the equations 

(17) and (18) . Thus, the matching computed at a coarser level 
of resolution is used in subsequent calculation of a matching. 
This is called a vertical reference between different levels. 
If m=0, there is no coarser level and this exceptional case 

15 will be described using Fig. 13. 

A horizontal reference within the same level is also 
performed. As indicated by the equation (20) in [1.3.3], 
£(m,3)^ £(m,2) f(m,i) respectively determined so as to be 

analogous to ff"''^), f^'D and f ("''O' . This is because a 

20 situation in which the submappings are totally different seems 
unnatural even though the type of critical points differs so 
long as the critical points are originally included in the 
same source and destination images. As can been seen from the 
equation (20), the closer the submappings are to each other. 



51 

the smaller the • energy becomes, so that the matching is then 
considered more satisfactory. 

As for f^™'°', which is to be initially determined, a 
coarser level by one may be referred to since there is no 
5 other submapping at the same level to be referred to as shown 
in the equation (19) . In this base technology, however, a 
procedure is adopted such that after the submappings were 
obtained up to f'"''^', f(m,o) recalculated once utilizing the 
thus obtained subamppings as a constraint. This procedure is 

10 equivalent to a process in which s=4 is substituted into the 
equation (20) and f^'"'^' is set to f'""'"^ anew. The above 
process is employed to avoid the tendency in which the degree 
of association between f*"''^' and f'"''^' becomes too low. This 
scheme actually produced a preferable result. In addition to 

15 this scheme, the submappings are shuffled in the experiment as 
described in [1.7,1], so as to closely maintain the degrees of 
association among submappings which are originally determined 
independently for each type of critical point. Furthermore, 
in order to prevent the tendency of being dependent on the 

20 starting point in the process, the location thereof is changed 
according to the value of s as described in [1.7]. 

Fig. 13 illustrates how the submapping is determined at 
the 0-th level. Since at the 0-th level each sub-image is 
consitituted by a single pixel, the four submappings f<°'^' are 



52 

automatically chosen as the identity mapping. Fig. 14 shows 
how the submappings are determined at the first level. At the 
first level, each of the sub-images is constituted of four 
pixels, which are indicated by solid lines. When a 
5 corresponding point (pixel) of the point (pixel) x in p'^'^* is 
searched within q*-*"^', the following procedure is adopted: 
1. An upper left point a, an upper right point b, a lower left 
point c and a lower right point d with respect to the point x 
are obtained at the first level of resolution. 

10 2. Pixels to which the points a to d belong at a coarser level 
by one, i.e., the 0-th level, are searched. In Fig. 14, the 
points a to d belong to the pixels A to D, respectively. 
However, the pixels A to C are virtual pixels which do not 
exist in reality. 

15 3. The corresponding points A' to D' of the pixels A to D, 

which have already been defined at the 0-th level, are plotted 
in q*-"-'^' . The pixels A' to C are virtual pixels and regarded 
to be located at the same positions as the pixels A to C, 
4. The corresponding point a' to the point a in the pixel A is 

20 regarded as being located inside the pixel A' , and the point 

a' is plotted. Then, it is assumed that the position occupied 
by the point a in the pixel A (in this case, positioned at the 
lower right) is the same as the position occupied by the point 
a' in the pixel A' . 



5. The corresponding points b' to d' are plotted by using the 
same method as the above 4 so as to produce an inherited 
quadrilateral defined by the points a' to d' . 

6. The corresponding point x' of the point x is searched such 
5 that the energy becomes minimum in the inherited 

quadrilateral. Candidate corresponding points x' may be 
limited to the pixels, for instance, whose centers are 
included in the inherited quadrilateral. In the case shown in 

Q. Fig. 14, the four pixels all become candidates. 

m 

jiy 10 The above described is a procedure for determining the 

m 

O corresponding point of a given point x. The same processing 

m 

" is performed on all other points so as to determine the 

submappings. As the inherited quadrilateral is expected to 
2 become deformed at the upper levels (higher than the second 



m 



level), the pixels A' to D' will be positioned apart from one 
another as shown in Fig. 3. 

Once the four submappings at the m-th level are 
determined in this manner, m is incremented (S22 in Fig. 12) . 
Then, when it is confirmed that m does not exceed n (323) , 
return to S21. Thereafter, every time the process returns to 
S21, submappings at a finer level of resolution are obtained 
until the process finally returns to S21 at which time the 
mapping f'"' at the n-th level is determined. This mapping is 
denoted as f"^'(r7=0) because it has been determined relative 
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to r}=0. 

Next, to obtain the mapping with respect to other 
different r] , ?? is shifted by Arj and m is reset to zero 
(S24) . After confirming that new r? does not exceed a 
5 predetermined search-stop value 77max(S25), the process returns 
to S21 and the mapping f'"' (tj = Arj) relative to the new r) is 
obtained. This process is repeated while obtaining 
f'''^\T] = ihTj) (z = 0,l,...) at S21. When r? exceeds rj^ax/ the process 
proceeds to S26 and the optimal 77 = »7opt is determined using a 
10 method described later, so as to let f'"M?7 = ^opt) be the final 
mapping f . 

Fig. 15 is a flowchart showing the details of the 
process of S21 shown in Fig. 12. According to this flowchart, 
the submappings at the m-th level are determined for a certain 
15 predetermined rj - In this base technology, when determining 

the mappings, the optimal A is defined independently for each 
submapping. 

Referring to Fig. 15, s and A are first reset to zero 

(5210) . Then, obtained is the submapping f = ' that minimizes 
20 the energy with respect to the then A (and, implicitly, r? ) 

(5211) , and the thus obtained submapping is denoted as 
£(in,s) ^ ^ ^ In order to obtain the mapping with respect to 
other different A, A is shifted by AX. After confirming 
that the new A does not exceed a predetermined search-stop 



value Aniax (S213), the process returns to S211 and the mapping 
£(m,s) (;j, = A>i) relative to the new A is obtained. This process 
is repeated while obtaining f (;i = /AA) 0" = 0,1,—) • When A 
exceeds X^^, the process proceeds to S214 and the optimal 
5 A = Aopt is determined , so as to let f^"'(A = Aopt) be the final 
mapping f (3214) . 

Next, in order to obtain other submappings at the same 
level, A is reset to zero and s is incremented (S215) . After 
confirming that s does not exceed 4 (S216) , return to S211. 

10 When s=4, f<"''°^ is renewed utilizing f*'"'^' as described above 
and a submapping at that level is determined. 

Fig. 16 shows the behavior of the energy Cj'"'*^ 
corresponding to f (>i = rAA) (/ = 0,1,...) for a certain m and s 
while varying A. As described in [1.4], as A increases, 

15 C^/''^ normally decreases but changes to increase after A 

exceeds the optimal value. In this base technology, A in 
which Cj™'*^ becomes the minima is defined as A opt- As observed 
in Fig. 16, even if Cf''"^ begins to decrease again in the range 
A>Aopt / the mapping will not be as good. For this reason, it 

20 suffices to pay attention to the first occurring minima value. 
In this base technology, A opt is independently determined for 
each submapping including f ''^^ . 

Fig. 17 shows the behavior of the energy Cf^ 
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corresponding to f'-^^(j] = iAr])(i = 0,\,—) while varying rj . Here 
too, Cj"^ normally decreases as ii increases, but Cj"^ changes 
to increase after rj exceeds the optimal value. Thus, r? in 
which Cj"^ becomes the minima is defined as 77 opt- Fig. 17 can 
5 be considered as an enlarged graph around zero along the 

horizontal axis shown in Fig, 4. Once 77 opt is determined, f^'^' 
can be finally determined. 

As described above, this base technology provides 
various merits. First, since there is no need to detect 

10 edges, problems in connection with the conventional techniques 
of the edge detection type are solved. Furthermore, prior 
knowledge about objects included in an image is not 
necessitated, thus automatic detection of corresponding points 
is achieved. Using the critical point filter, it is possible 

15 to preserve intensity and locations of critical points even at 
a coarse level of resolution, thus being extremely 
advantageous when applied to object recognition, 
characteristic extraction, and image matching. As a result, 
it is possible to construct an image processing system which 

20 significantly reduces manual labor. 

Some further extensions to or modifications of the 
above-described base technology may be made as follows: 
(1) Parameters are automatically determined when the matching 
is computed between the source and destination hierarchical 
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images in the base technology. This method can be applied not 
only to the calculation of the matching between the 
hierarchical images but also to computing the matching between 
two images in general . 
5 For instance, an energy Eq relative to a difference in 

the intensity of pixels and an energy Ei relative to a 
positional displacement of pixels between two images may be 
used as evaluation equations, and a linear sum of these 
equations, i.e., Etot=aEo+Ei, may be used as a combined 
10 evaluation equation. While paying attention to the 

neighborhood of the extrema in this combined evaluation 
equation, a is automatically determined. Namely, mappings 
which minimize Etot are obtained for various a's. Among such 
mappings, a at which Etot takes the minimum value is defined 
15 as an optimal parameter. The mapping corresponding to this 
parameter is finally regarded as the optimal mapping between 
the two images. 

Many other methods are available in the course of 
setting up evaluation equations. For instance, a term which 
20 becomes larger as the evaluation result becomes more 

favorable, such as 1/Ei and I/E2, may be employed. A combined 
evaluation equation is not necessarily a linear sum, but an n- 
powered sum (n=2, 1/2, -1, -2, etc.), a polynomial or an 
arbitrary function may be employed when appropriate. 
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The system may employ a single parameter such as the 
above a , two parameters such as r) and A as in the base 
technology, or more than two parameters. When there are more 
than three parameters used, they may be determined while 
5. changing one at a time. 

(2) In the base technology, a parameter is determined in a 
two-step process. That is, in such a manner that a point at 
which Cf''^ takes the minima is detected after a mapping such 
that the value of the combined evaluation equation becomes 

10 minimum is determined. However, instead of this two-step 

processing, a parameter may be effectively determined, as the 
case may be, in a manner such that the minimum value of a 
combined evaluation equation becomes minimum. In this case, 
aEo+;8Ei, for example, may be used as the combined evaluation 

15 equation, where a + i8=l may be imposed as a constraint so as 
to equally treat each evaluation equation. The automatic 
determination of a parameter is effective when determining the 
parameter such that the energy becomes minimum. 

(3) In the base technology, four types of submappings 

20 related to four types of critical points are generated at each 
level of resolution. However, one, two, or three types among 
the four types may be selectively used. For instance, if 
there exists only one bright point in an image, generation of 
hierarchical images based solely on f'"^'^' related to a maxima 
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point can be effective to a certain degree. In this case, no 
other submapping is necessary at the same level, thus the 
amount of computation relative on s is effectively reduced. 
(4) In the base technology, as the level of resolution of an 
5 image advances by one through a critical point filter, the 
number of pixels becomes 1/4. However, it is possible to 
suppose that one block consists of 3X3 pixels and critical 
points are searched in this 3X3 block, then the niomber of 
pixels will be 1/9 as the level advances by one. 

10 (5) In the base technology, if the source and the 

destination images are color images, they would generally 
first be converted to monochrome images, and the mappings then 
computed. The source color images may then be transformed by 
using the mappings thus obtained. However, as an alternate 

15 method, the submappings may be computed regarding each RGB 
component . 

Preferred Embodiments for Image Data Coding and Decoding 

In the above-described base technology, the 
20 correspondence data are generated by computing a matching 
between key frames and, based on this correspondence 
information, an intermediate frame is generated. As such, 
this technology can be used for the compression of motion or 
moving pictures. In fact, experiments are beginning to show 
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evidence of both picture quality and compression rates that 
are superior to those of MPEG. An image coding technology 
utilizing the base technology will now be described with 
reference to Figs 18 to 34. This image coding technology, 
which takes transmission errors into consideration, is 
intended to be a useful elemental technique for the 
distribution of motion pictures utilizing the base technology. 

(Coding Side) 

Fig. 18 shows a structure of an image coding apparatus 
10 according to an embodiment of the invention. The image 
coding apparatus 10 includes an image input unit 12 which 
receives data of key frames KFi, a matching processor 14 which 
generates a correspondence data file or files Ci,j by computing 
a matching between key frames KFi and KFj, and a stream 
generator 16 which generates a coded data stream CBS (Coded 
Bit Stream) (sometimes simply referred to as a data stream 
hereinafter) , by incorporating the key frame data and the 
correspondence data file or files. The data stream may be 
stored in a storage unit 18 as needed. 

Key frames are defined, for instance, by extracting 
frames of motion pictures at predetermined time intervals, 
such as 0,5 seconds, and intermediate frames, which are the 
frames between key frames, are generated by an interpolation 
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computation based on the correspondence data file at a 
decoding stage. It is to be appreciated here that the image 
input unit 12 may either receive already existing key frames 
from an external storage device (not shown) , network (not 
shown) or the like or may be an image photographing device 
which captures images, such as a digital camera. 

The matching processor 14 performs a pixel-by-pixel 
matching between two key frames, based on critical points or 
otherwise, applying the base technology or some other 
arbitrary technology. Unlike the case with the base 
technology, the matching processor 14 computes a matching not 
only between adjacent key frames but also key frames that are 
a predetermined distance away from each other (i.e. 
predetermined at timings) . This matching will be called a 
"spare matching" hereinafter. The spare matching generates 
correspondence data for coping with transmission errors or the 
like. Namely, when there is a data error in key frame data or 
in a correspondence data file during transmission and so 
forth, this will be remedied by the spare matching 
correspondence data file (also referred to as a "'spare file" 
hereinafter) which is obtained by the spare matching. 

Now, for nonnegative integers i and j, let us define 
data of an i-th key frame and a j-th key frame as KFi and KFj, 
respectively, and define a correspondence data file between 
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the i-th key frame and the j-th key frame as Ci,j. Then, the 
matching processor 14 can be configured in a manner such that 
the matching processor 14 generates Ci, i+i by utilizing KFi and 
KFi+i among input key frame data, and generates a spare Ci,j by 
5 utilizing KFi and KFj {j>i + l) . It is to be noted that Ci,j is 
generated with lower frequency compared to Ci, i+i. In order to 
check for transmission errors, the matching processor 14 adds 
a parity bit to the correspondence data files Ci,i+i and Ci,j. 
In order to check for transmission errors in the key frames, 

10 the image input unit 12 or the stream generator 16 may add a 
parity bit or bits to the key frames, as described in more 
detail below. 

Fig. 19 shows a relationship between the ordinary 
correspondence data file (also referred to as a "main file" 

15 hereinafter) and a spare file. Referring to Fig. 19, the main 
file is generated one at a time between adjacent key frames in 
sequence. On the other hand, spare files Co,6, Ce,i2 - are 
generated with some interval therebetween. 

Fig. 20 shows a structure of the data stream CBS which 

20 is generated in a manner that the main files and spare files 
shown in Fig. 19 are incorporated thereto by the stream 
generator 16. If no countermeasure against error is 
considered, this stream will be such that key frames and main 
files appear in a predetermined sequence. Specifically, the 
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key frames and the main files appear in the order of KFq, KFi, 
Co,ir KF2, Ci,2, KF3, C2,3/-. The spare files Co,6, C6,i2 are 
inserted before KFo and KF7, respectively. A method for 
coping with transmission error using the spare file will be 
described below when describing a decoding stage. 

Although the insertion position of the spare file is 
somewhat arbitrary, the following two points need to be taken 
into consideration in general. As a first point, it is 
preferable that key frames and main files covered by a spare 
file (also simply referred to as a ^'covering area" 
hereinafter) are, to a certain degree, away from the position 
of the spare file. This is because a case when the 
transmission error extends or occurs simultaneously to both 
the spare file and the cover file should be avoided. As a 
second point, though contrary to the first point, the both 
shall be disposed close to each other, to a certain degree. 
This is because the spare file needs to be present at the time 
of reproduction of the covering area, and the spare files 
cannot be sent later unless the covering area is in a stand-by 
state. On the other hand, if the spare files are sent in 
advance, these need to be stored in a buffer for a long period 
of time. Thus, it is desirable that a capacity of a device 
and so forth be taken into account, and the positions of the 
spare files be determined through some experiments or the 
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like. It is to be noted that the same spare file may be 
embedded a plurality of times in different places. 

Fig- 21 shows another possible relationship between a 
main file and a spare file. Here, Co, 6/ €4,10, etc. are 
5 generated in a manner that the covering areas of the spare 
files are overlapping with one another. Fig. 22 shows a 
structure of the data stream CBS which corresponds to Fig. 21. 
Here, spare files €4,10 and C8,i4/ ... are inserted before KF5, 
KF9, respectively. The error tolerance further improves 

10 due to the overlapping effects. 

Fig. 23 shows still another relationship between a main 
file and a spare file. Here, in addition to the spare files 
Co,6r C.s,i2,—r Spare files (Co, 12/- Ci2,24f ■■• ) which cover a longer 
interval are added. Fig. 24 shows a structure of the data 

15 stream CBS which corresponds to Fig. 23. Here, spare files 
Co, 12/ Co,6f C6,i2/- are inserted in the head of the stream and 
before KF3 and KF7, respectively. In a case of this 
structure, when there is an error in, for example, the 
covering area of Co, 6/ Co, e will be used at a primary stage. If 

20 there is also an error in Co, 6/ then Co, 12 can be used at a 

secondary stage, thus further improving the error tolerance. 

As described above, error checking data, such as parity 
bits, are added to key frames and correspondence data files. 
For key frames, it may not be necessary to check for errors in 
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the whole frame since some types of errors, such as those at 
the edges, may be less important. Fig. 25 shows an area of a 
key frame KFi in which a data error is to be detected. Here, 
sixteen block regions Ro - R15 are provided in the vicinity of 
the center of an image, and one-dimensional or two-dimensional 
parity bit is added for each block. In this case, it has been 
determined that the area close to the center of the key frame 
is of higher importance and only errors in that area will be 
dealt with. However, of course, the whole area of the image 
may also be taken into account. Moreover, a region considered 
to have high importance in general such as a region having a 
high image intensity may be dynamically determined, so that 
errors in said region may be considered. 

Fig. 2 6 shows conceptually an area of a correspondence 
data file Ci,j in which a data error is to be detected. This 
file generally stores text data and, here, is simply divided 
into a plurality of data blocks To - T4, so that each block is 
given a parity bit. 

(Decoding Side) 

Fig. 27 shows a structure of an image decoding apparatus 
40 according to an embodiment. The image decoding apparatus 
40 includes a stream input unit 42 which receives or acquires 
a coded data stream CBS, an intermediate image generator 44 



which generates an intermediate frame from the stream by 
interpolation based on, for example, the base technology, and 
a display controller 50 which performs a processing to 
display key frames and intermediate frames as moving pictures. 
Display data generated by the display controller 50 are 
outputted to a display device where images are reproduced. An 
error detector 4 6 detects an error in the data stream. If 
there is an error in the data stream, the error detector 46 
notifies an error controller 48 to that effect. The error 
controller 48 incorporates a predetermined error avoidance 
processing to generation of intermediate frames in the 
intermediate image generator 44. The stream input unit 42 and 
the intermediate image generator 4 4 have respective buffer 
memories or a shared buffer memory which are or is utilized as 
a work area for data necessary for decoding or error avoidance 
processing. Hereinafter, data of an intermediate frame for 
KFi and KFj will be denoted as IFi,j. 

Fig. 28 illustrates a principle of generation of 
intermediate frames by the intermediate image generator 44. 
This processing itself was also introduced in the base 
technology. Namely, when there are two frames KF5 and KFe and 
their main file 05,5, it is first confirmed by this main file 
C5,6 that a point po (xo, yo) on the key frame KF5 corresponds to 
a point pi(xi, yi) on other key frame KFe. A point is 



67 

specified by a coordinate (xi, yi) and a pixel value pi. Thus, 
a point pt(xt, yt) on an intermediate frame IF5,6(t) for these 
key frames is obtained by interpolating the coordinate and the 
pixel value on the time axis. 

Fig. 29 shows an error avoidance processing controlled 
by the error controller 48 when an error is detected in KFe of 
a data stream CBS that is being processed by the image 
decoding apparatus 4 0'. Here, a key frame KF5, which is a key 
frame immediately before KFg, is copied to KFg, and this 
copied KF5 is utilized as KFs. Of course, alternatively a key 
frame which comes after KFe may be utilized, or a weighted 
average or the like of key frames before and after KFg may be 
utilized. 

Fig. 30 shows another error avoidance processing when an 
error is similarly detected in KFg. Here, the error- 
containing KFe is not used at all. Instead, each point of KF5 
is moved using €5,5, so as to generate IF5,6. In other words, 
movement of pixels only is traced by 05,6/ and the change of 
pixel values to be traced primarily is omitted because KFg is 
not available. IF5,6 which is the intermediate frame 
representing part of KFg is generated by a transformation or 
deformation of KF5. 

Fig. 31 shows an error avoidance processing when an 
error is detected in two key frames KFe and KF7, in a row. 



Here, KFe- and KF7 are not used, instead, each point of KF5 is 
moved by utilizing 05,5, so as to generate IF5,6. As a result, 
KFe is virtually generated as a destination image. 
Thereafter, each point of this virtually generated KFe is 
moved by utilizing Cg,?, so as to generate IFg,,. Then, KF7 is 
virtually generated as a destination image, and is utilized 
for subsequent processing. 

Fig. 32 shows still another type of error avoidance 
processing. Here, an error is detected in more than one of 
key frames KFi - KF5, or, as the case may be, a main file 
related thereto, for example, 02,3- Based on an instruction 
from the error controller 48, the intermediate image generator 
44 first selects a spare file generated for key frames that 
interpose an error portion, among spare files included in the 
data stream CBS. Referring to Fig. 32, Co,6 is detected, thus 
this Co, 6/ KFo and KFe are utilized. As a result, intermediate 
frames IFo,6 over a rather long interval therebetween is 
generated. Even if the error continues over a relatively long 
period of time, this method of reproducing images with a 
certain level of image quality is extremely useful. 

The following describes a method which generalizes the 
above example error avoidance processes. Namely, when an 
error occurs in KFi+i, an error avoidance processing is 
performed in a manner such that Ca,b, where a ^ i+1, 1+2 < b or 



a < i+1, i+2 S b, is detected among spare files, and IFa,b is 
generated by utilizing Ca,br KFa and KFb, thereby this IFa,b or a 
part thereof is substituted for IFi+i,i+2. 

There may also be a case where, Ca,b, a < i+1, i+2 < b, 
is used. This may occur when such a spare file is the only 
one that is available or when the number of spare files is to 
be reduced deliberately. 

Fig. 33 shows still another error avoidance processing. 
Here, suppose that error occurred in a main file Cb,^. In this 
case, C4,5 which is the main file immediately before Cs^g is 
copied and is substituted for the error-containing Cs^q. Of 
course, alternatively, a main file which comes after C5,6 may 
be copied, or a weighted average or the like of main files 
before and after C5,6 may be utilized. 

Fig. 34 shows still another error avoidance processing. 
Here, too, an error occures in the main file C5,6. Here, the 
content of the file 05,5 is replaced with correspondence data 
that indicates an identity mapping. The "identity mapping" is 
a mapping such that f (x) = x. Here, the identity mapping is 
thought of as a case where all of the pixels of KF5 correspond 
to the same coordinates of KFe- This method is effective 
because, in a movie or video sequence of frames, differences 
between two adjacent key frames are, in general, not very 
large. This method may also be effective in a case where main 
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files are lost in a comparatively consecutive manner. 

It is useful to note that, even when error is detected 
in a data stream CBS, the execution of the error avoidance 
processing may be abandoned, put off or reconsidered if it is 
5 determined that the seriousness or importance of the error is 
substantially low. Recall that, as shown in Fig. 25, regions 
of a key frame may be configured so that error checking may be 
M performed primarily on the important area or areas only, 

however, even though error checking is possible over the whole 
10 area, errors may be ignored at the decoding side if the error 

^* occurs in, for example, an edge of an image. In particular, 

e 

it may be possible for some cooperation and special 
W arrangement to be effected between the coding side and 

J' decoding side to determine which types of errors are to be 

'"*15 ignored or corrected. 

Although the present invention has been described by way 
of exemplary embodiments, it should be understood that many 
changes and substitutions may be made by those skilled in the 
art without departing from the scope of the present invention 
20 which is defined by the appended claims. Thus, the present 

invention is not limited to the embodiments described herein, 
and various modifications thereto are also effective as and 
encompassed by the present invention. 



